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Abstract 

Multidimensional minimax single facility location problems with 
Chebyshev distance are examined within the framework of idempotent 
algebra. A new algebraic solution based on an extremal property of 
the eigenvalues of irreducible matrices is given. The solution reduces 
both unconstrained and constrained location problems to evaluation 
of the eigenvalue and eigenvectors of an appropriate matrix. 
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1 Introduction 

Location problems form one of the classical research domains in optimiza- 
tion that has its origin dating back to XVIIth century. Over many years 
a large body of research on this topic contributed to the development in 
various areas including integer programming, combinatorial and graph opti- 
mization. Among other solution approaches to location problems are models 
and methods of idempotent algebra [U O El [H El E] , which find expanding ap- 
plications in the analysis of actual problems in engineering, manufacturing, 
information technology, and other fields. Specifically, an algebraic solution 
to a one-dimensional location problem on a graph is proposed in |2]. A 
constrained location problem and its representation in terms of idempotent 
algebra are examined in [8J [9] . 

In this paper, we consider a multidimensional minimax single facility 
location problem with Chebyshev distance, and show how the problem can 
be solved based on new results of the spectral theory in idempotent algebra. 
The aim of the paper is twofold: first, to give a new algebraic solution 
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to the location problem, and second, to extend the area of application of 
idempotent algebra. The rest of the paper is as follows. We begin with 
an overview of preliminary algebraic definitions and results. Specifically, 
we present an extremal property of the eigenvalue of irreducible matrices. 
Furthermore, we examine an unconstrained minimax location problem and 
represent it in the terms of idempotent algebra. A new solution is given that 
reduces the problem to evaluation of the eigenvalue and eigenvectors of an 
irreducible matrix. Finally, the solution is extended to solve a constrained 
location problem. 

2 Preliminary Results 

We start with a brief overview of definitions, notation and preliminary re- 
sults of idempotent algebra that underlie the solution approach developed 
in subsequent sections. Further details can be found in [H El El IH El [6] . 

2.1 Idempotent Semifield 

Let X be a set with two operations, addition © and multiplication ® , and 
their respective neutral elements, and 1. We suppose that (X, 0, 1, ffi, ®) 
is a commutative semiring where the addition is idempotent and the mul- 
tiplication is invertible. Since the nonzero elements of the semiring form a 
group under multiplication, the semiring is usually referred to as idempotent 
semifield. 

The integer power is defined in the ordinary way. Let us put X + = 
X \ {0}. For any x 6 X + and integer p > 0, we have x° = 1, x p = 
xP- 1 ®x = x®x p ~ 1 , x~P = (x~ l Y, and P = 0. 

We assume that the integer power can naturally be extended to the case 
of rational and real exponents. 

In what follows, we drop, as is customary, the multiplication sign ®. 
The power notation is used in the sense of idempotent algebra. 

The idempotent addition allows one to define a relation of partial order 
< such that x < y if and only if x © y = y . From the definition it follows 
that x < x(By and y < well as that the addition and multiplication 

are both isotonic. Below the relation symbols are thought of as referring to 
this partial order. 

It is easy to verify that the binomial identity now takes the form (x © 
y ) a = x a ®y a for all a > 0. 

As an example, one can consider the idempotent semifield of real num- 
bers 

Kmax,+ = (RU {-oo}, -oo, 0, max, +). 

In the semifield R maXj+ , there are the null and identity elements defined 
as = — oo and 1 = 0. For each i£R, there exists its inverse x _1 equal 
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to the opposite number — x in conventional arithmetic. For any x,y G R, 
the power x y coincides with the arithmetic product xy . 

2.2 Vectors and Matrices 

Vector and matrix operations are routinely introduced based on the scalar 
addition and multiplication defined on X . Consider the Cartesian power X ra 
with its elements represented as column vectors. For any two vectors x = 
(xj) and y = (yi) , and a scalar c G X, vector addition and multiplication 
by scalars follow the rules 

{x © y\i = Xi © yi, {c£c}j = cxj. 

The set X ra with these operations is a vector semimodule over the idem- 
potent semifield X. 

As it usually is, a vector y G X ra is linearly dependent on vectors 
x ii • • • j x m G X n , if there are scalars ci, . . . , c m G X such that y = ci»i © 
• • • © c m x m . In particular, the vector y is collinear with a:, if y = ex. 

For any column vector x = (xj) G X™ , we define a row vector x = 
(x~) with its elements x~ = x~ x . For all x,y G X™ , the componentwise 
inequality x < y implies x > y . 

For any conforming matrices A = (ay), B = and C = (cy), 

matrix addition and multiplication together with multiplication by a scalar 
c G X are performed according to the formulas 

{^4 © B}ij = ciij © bij, {BC}ij = ^^^ifcCfcj, {c^4}jj = ca^-. 

A matrix with all zero entries is referred to as zero matrix and denoted 
by 0. 

Consider the set of square matrices X nxn . The matrix that has all 
diagonal entries equal to 1 and off-diagonal entries equal to is the identity 
matrix denoted by /. 

With respect to matrix addition and multiplication, the set X nxn forms 
idempotent semiring with identity. 

For any matrix A ^ and an integer p > 0, we have A = I and 
A p = ap- 1 A = AA 13 - 1 . 

The trace of the matrix A = (ay) is defined as tr A = an © • • • © a nn . 

A matrix is irreducible if and only if it cannot be put in a block triangular 
form by simultaneous permutations of rows and columns. Otherwise the 
matrix is reducible. 

2.3 Eigenvalues and Eigenvectors of Matrices 

A scalar A is an eigenvalue of a square matrix A G X raxra if there exists a 
vector x G X n \ {0} such that Ax = Xx. Any vector x ^ that satisfies 
the equation is an eigenvector of A , corresponding to A . 
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If the matrix A is irreducible, then it has only one eigenvalue given by 

n 

A = t?/ m {A m ). (1) 

m=l 

The corresponding eigenvectors of A have no zero entries and are found 
as follows. First we evaluate the matrix 

a x =A- 1 Ae-e(A" 1 A) n 

Let af and be respective column i and diagonal entry (i,i) of A x . 
Consider the subset of columns a* such that = 1, i = 1, . . . , n. In the 
subset, find those columns that are linearly independent of the others, and 
take them to form a matrix A + . 

The set of all eigenvectors of A corresponding to A (together with zero 
vector) coincides with the linear span of the columns of A x , whereas each 
vector takes the form 

x = A + v, 

where v is a nonzero vector of appropriate size. 

3 An Extremal Property of Eigenvalues 

Suppose A 6 X nxn is an irreducible matrix with an eigenvalue A. Consider 
a function (p(x) = x^ Ax defined on X™ . It has been shown in |1U^ 111] that 
(p{x) has A as its minimum, which is attained at any eigenvector of A. 

Now we improve the above result by extending the set of vectors that 
provide for the minimum of (p{x) . 

Lemma 1. Let A = (ciij) 6 x nxn be an irreducible matrix with an eigen- 
value A. Suppose u = (ui) and v = {yi) are eigenvectors of the respective 
matrices A and A T . Then it holds that 

min x~ Ax = A, 

where the minimum is attained at any vector x = (u^v . . . , u^v^~ 1 ) T 
for all a such that < a < 1 . 

Proof. It is easy to verify as in [11] that any vector x with nonzero elements 
satisfies the inequality x~ Ax > A. Indeed, since xu > (a^w) -1 / , we have 
x Ax = x Axu u > x Au(x~ u)~ x = Xx^u(x^u)^ 1 = A. 

It remains to present a vector x that turns the inequality into an equality. 
With x = u we immediately have x~Ax = \u~u = A. Similarly, if 
x = (v~) T , then x~ Ax = x T A T (x~) T = v~A T v = \v~v = A. 
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Let us now take any vector x = (xi) with elements Xi = ufv^ 1 , where 
a is a real number such that < a < 1 . With the following calculations 

n n n n 

A = {u-AuY{v-A T vf-« = ©©V a ^©0<" 1 «" a 

i=l j=l k=l 1=1 

n n n n 

i=i j=i i=i j=i 

we arrive at the inequality x Ax < A. Since the opposite inequality is 
always valid, we conclude that the vector x satisfies the condition x^ Ax = 
A. □ 



4 The Unconstrained Location Problem 

In this section we consider a minimax single facility location problem with 
Chebyshev distance. For any two vectors r = (n,...,r n ) T and s = 
(si, . . . ,s n ) T in R n , the Chebyshev distance (Loo ° r maximum metric) is 
calculated as 

p(r,s)= max |r» — Sj|. (2) 

l<i<n 

Given m > 2 vectors r, = (rij, . . . , r n i) T € R n and constants Wi £ R, 
z = 1, ...,m, the problem under consideration is to find a vector x = 
(xi, . . . , x n ) T so as to provide 

min max cc) + wA (3) 

u€R" l<i<m 

It is not difficult to solve the problem by using geometric arguments 
(see, eg, [12], Q3]). Below we give a new algebraic solution that is based on 
representation of the problem in terms of the idempotent semifield R m ax,+ 
and application of the result from the previous section. First we rewrite ([2]) 
as follows 

p(r, s) = s r © r~s. 
Denote the objective function of the problem by <p(x) . With the vectors 

p = w\ri © • • • © w m r m , q = wir^ © • • • © w m r m , 

we can write 

m m 

fi x ) = @WiP(ri,x) = (Qwi(x~ri © r^x) = x p®q x, 

i=l i=l 

and then represent problem ([3]) as 

min ip(x) = min (x~p © q~x). (4) 
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Furthermore, we introduce a vector y and a matrix A of order n + 1 as 
follows 

-(:)■ Hi q « 

Since we now have </?(a:) = a? p@q x = y Ay , problem ([5]) reduces to 
that of the form 

min y Ay. (5) 

Note that the vectors y £ R n+1 that solve © do not always have an 
appropriate form to give a solutions to (jU). Specifically, to be consistent to 
d3]) , the vector y must have the first element equal to 1 . 



5 Algebraic Solution of the Unconstrained Prob- 
lem 

Consider problem ([5]), and note that the matrix A is irreducible. It follows 
from Lemma [1] that the minimum in ([5]) is equal to the eigenvalue A of the 
matrix A. For all k = 1, 2, . . . we have 

and therefore, ti{A 2k ~ l ) = 0, tr(A 2fc ) = (q~p) k ■ Finally, application of (P) 
gives 

A = (q-p) 1/2 . 

To find vectors that produce the minimum in (0), we need to derive the 
eigenvectors of the matrices A and A T . Note that A T is obtained from A 
by replacement of p with (q~) T and q~ with p T . Therefore, it will suffice 
to find the eigenvectors for A , and then turn them into the eigenvectors for 
A T by the above replacement. 

To get the eigenvectors of A, we first find the matrix A x . Since for any 
k = 1, 2, . . . it holds that 

(A-'A)*- 1 = (,-p)-'/ 2 ( * «" 

<*-'A)» = (.-pr'(Y 

we arrive at the matrix A x in the form 

a* = x-'a e • • • e (A-^r +1 = ( , f 1/2 ^l p \'_T Q - 

V (9 p) 7 p (q p) pq 

It is not difficult to verify that in the matrix A x , any column that has 1 
on the diagonal is collinear with the first column. Indeed, suppose that the 
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submatrix (q~p)~ PQ~ has a diagonal element equal to 1, say the element 
in its first column (that corresponds to the second column of A x ). In this 
case, we have q~p = q : [ 1 pi , whereas the matrix A x takes the form 



A" 



t 



1/2 -1/2 _ 

q{ p 1 q 

1/2 -1/2 -1 

Qt Pi p qwi pq 



1/2 -1/2 

q{ p 1 p 



-1/2 -1/2 

Qi Pi 
Pi X P 



where the second column proves to be collinear with the first. 

Let us construct a matrix A + that includes such columns of A x that 
have the diagonal element equal to 1. Since all these columns are collinear 
with the first one, they can be omitted. With the matrix A + formed from 
the first column of A x , we finally represent any eigenvector of A as 

1 

(q-p)- l/2 p 

By replacing p with (q~) T and q with p T , we get the eigenvectors of 



u 



s e 



A T 



1 

{q-p)- 1/2 {q-) T 



t, 



t e 



It follows from Lemma Q] that the solution of (j3|) takes the form 
\ / 1 \ 



y 



^0 



/ u v i 

a„.a— 1 



1 u l 



V u n+l v n+l J 



(q-p)^ 2 - a p?qt a 



\ (q-p) l/2 - a P> l r r a 



s a t a ~\ 



s,te R. 



With the condition that the first element of y must be equal to 1, 
we have to set s a t a ~ l = 1. Going back to problem we arrive at the 
following result. 

Lemma 2. The minimum in problem @ is given by X = (q^p) 1 ^ 2 , and it 
is attained at any vector 

( P?ql- a 



x = A 



l-2a 



< a < 1. 



V P a n^n a 

With the usual notation, we can reformulate the statement of Lemma [2] 
as follows. 



Corollary 3. The minimum in ([3]) is given by A = max(j>i — qi, . 
q n )/2, and it is attained at any vector 



< a < 1, 



,Pn 





f Pi-X\ 




( ?i + A \ 


x = a 


^ Pn ~ A J 


+ (1 - a) 


V Qn + a / 



where pi = max(rji + wi, 
for each i = 1, . . . , n. 



, r im + Wm) , qi = min(rji - wi, 



W r . 
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6 A Constrained Location Problem 



Suppose that there is a set S G R n given to specify feasible location points. 
Consider a constrained location problem that is represented in terms of the 
semifield K m ax,+ in the form 

m 

min ff^Wip(ri,x). (6) 

To solve the problem we put it in the form of (jU) by including the area 
constraints into the objective function of an unconstrained problem. First we 
slightly transform problem to enable accommodation of the constraints 
in a natural way. With the notation 

Po = win © • • • © w m r m , q = win • • • w m r m , 

Ao = {q Po) 1/2 , <Po{x) = Aq 1 {x~Pq q x), 

we turn to the problem 

min lpq(x). 

33GR™ 

It follows from Lemma [2] that the last problem has its minimum equal 
to 1 = 0, whereas its solution set obviously coincides with that of flU) . 

Suppose that there are constraints on the maximum distance from the fa- 
cility location point to each given points that determine the feasible location 
set in the form S = {x G R n \p(ri, x) < di,i = 1, ... , rn}. 

For each i = 1, . . . , to, the inequality x~ri © nx = p(n, x) < d{ can 
be rewritten in an equivalent form as d^ 1 x~ri © d^ 1 r^x < 1, whereas all 
constraints can be replaced with one inequality 

X (d-^ 7*i © • • • © d m ^rn) © (^i Ti © * * * © d m T rn ^X 1. 

We introduce the notation 

Pi = d^ V\ © ■ ■ ■ © d m T m , q^ = d-y v-y © • • • © d m v m: 
(fl(x) = x pi © q^x, 

and note that (pi(x) < 1 when and only when the above constraints are 
satisfied. 

Furthermore, we put 

p = Aq 1 po © pi , q~ = Aq 1 Qq © Qi , 

ip(x) = tpo(x) © = x~p © q x. 

Now we can replace the original problem ([6]) by an unconstrained prob- 
lem with the objective function <p(x) that has the form of problem (HJ with 
the solution given by Lemma [2j It is clear that both problems give the 
same solution set provided that a solution exists. At the same time, the 
new problem allows one to get approximate solutions in the case when the 
original problem does not have a solution. 
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